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ABSTRAC'I 

A ebulliometric techniaue has been "used 

to measure the vapor pres sure -hoi ling point data of Benzyl 
Benzoate, Brornoform, 1-Bromonaphthalene, Cycloheptane, 

Nonyl Alcohol (n-llonanol), Nonyl Phenol, Propylene Trimer 
from about go to 900 mm Hg. pressure. Differential 
ehulliometers which enable the measurement of boiling and 
condensing temperatures at a" particular pressure have been 
used. Y^ater has been used as the reference compound. 
Temperatures are measured to an accuracy of + 0,01°C 

and pressures to + 0.1 mm Hg, Literature data are available 
for Benzyl Benzoate, Brornoform, 1-Bromonaphthalene and Honyl 
Alcohol, but the data are old and unreliable. Actual 
accuracy of temperature and pressure for each compound .are 
reported. The data are fitted to the Antoine and Cox vapor 
pressure equations with good accuracy. 
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CHAPTBR 1 

BTTRODUCTIOR 

Out of a very large number of organic substances 
synthesized so far, reliable experimental vapor pressure- 
boiling point data are available for at most one thousand 
substances. Much of this published information is available 
in the Thermodynamic Research Center compilations ( 18 , 19)5 
Timmermans ( 24 ), Reid and Sherwood ( 16 ), Stull (20), Jordan ( 9)5 
and'Kanpur Critical Data Center compilation (IO). Experimental 
methods of measurement are elegantly discussed by Hala ,et,al. 
(6), In the present investigation, vapor pressure-boiling 
point measurements from, 90 to 900 mm mercury pressure were 
carried out with SY/ietoslawski type differential ebulliometers 
described in detail by Osborne and Douslin (12). This method 
was selected over the static methods because at the pressures 
under consideration here, differential ebulliometry gives very 
precise and accurate results. In the ebulliometric apparatus, 
measurements are made by directly comparing the boiling 
temperature of the sample with the boiling temperature of a 
suitable reference substance, pieferably water, in a common 
pressure manifold system’ blanketed with dry nitrogen. Tempera- 
tures were measured by copper-constantan thermocouples (3) 
calibrated against a Deeds and Rorthrup platinum resistance 
thermometer and K-3 potentiometer assembly. Pressures were 
roughly measured by a mercury manometer and accurately 
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calculated from 'botllng temperatures of the- referesce 
substance, water? the pressure value at each such temperature 
was calculated from the knovm vapor pres sure -boiling point 
data of vmter (18). The accuracy of the temperature and pressure 
measurements were expected to be + 0,02°C and 0,1 mm Hg, 
respectively. 

These vapor pressures are required in the prediction 
of vapor-liquid eqi;ilibria of multicomponent mixtures from 
the corresponding binaries. The compounds selected for the 
present investigation ?i;ere partly based upon our vapor-liquid 
equilibrium studies and partly because some high purity 
compounds were available. The following compounds were 
studied i 

1 . Benzyl benzoate' 

2 . Bromoform 

3. 1-Bro mo naphthalene 

4. Cycloheptane 

5. n-nonanol 

6. Nonyl phenol 

7. Propylene Trimer , 

Experimental data are available for benzyl benzoate (5), 

1-bro mo naphthalene (15), bromoform (l4)j and n-nonanol (20) 
but the data are old and unreliable. 

The experimental data were fitted to the Antoine (23) 
and the Cox (4) equations. The Antoine equation is 
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log^oP(mm) = 1 - ( 1 ) 

The totolne constants A,B,0 were determined by a least squares 
regression (l) with optimization method by Rosenbrock 
>-^s .1 (17) « The Antoine equation gives extremely good 

fit to the data in the range of 10 mm to 1500 mm Hg. pressure. 
The Cox equation 


log (P/760) = a (1 
cox 


log A = a + bT + cT^ 
■ cox 


T°K^ 


( 2 ) 


BP = normal boiling point; °K 

has been found to be better than the Antoine equation for 

fitting the data upto about 2000 mm Hg. pressure (4,12). The 

ideal Xi. values at the normal boiling point (at 760 mm Hg. 

pressure) were calculated by using Eqn. ( 1 ) and ■Eqn.(2) and 

the Clapeyron equation, Eqn, (5) assuming v v, . 
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dT 




V 
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.:MI P 

V 

RT"^ 


(3) 


The non-ideality correction could not be calculated because 
of the lack of necessary data. 
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CHAPTER 2 


E XPERIMEHTAE APPARATUS 

T1i 9 experluiental apparatus used in the present investi- 
gation is shcv’n in Pig. 1. The differential ehulliotneter 
shov/n in detail in Pig. 2, forms the heart of the apparatus. 

In this methods the boiling temperature of a substance is 
compared directly with that of water under nitrogen atmosphere. 

Pifferential Ebulliometer ; The ebulliometer used here is 
described by Osborne and Pouslin (12) and is a modification 
of the one designed by Swietoslawski (21), As shown in Pig. 2, 
it consists of a boiler and condenser sections provided with 
re-entrant wells for measurement of boiling and condensing 
temperatures. The boiler 16 is a glass tube of 3 cm. dia. 
and 11 cmd long provided with a glass re-entrant well of 
2 cm. dia for the heater. The liquid in the boiler is heated 
by means of an electric heater placed in the well containing 
silicone fluid (Iliiion Carbide) for a heat transmitting fluid. 
The heater is made of nichrome wire (24 gauge, 7,8 ohra/meter 
resistance) wound on a glass tube which"' is wrapped with 
asbestos clctlio Heat input to the heater is controlled through 
a variac. The heater is designed to cause necessary boiling 
rate corxesponding to the desired pressure of the liquid. 

To promote steady '-oiling, the inner surface of the well 
exposed to the sample was provided with glass beads. 
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Fig. 2 - Differential Ebulliometer, 
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The boiler is connected to the condenser portion 
with a percolation tube, 2^ of 1,3 cm. dia. and 16 cm, long. 

The ebulliometers are insulated to prevent heat loss, parti- 
cularly around the thermometer wells. However, the percolation 
tube remained essentially iminsulated (top 5 cm was insulated), 
so that any superheat in the boiling liquid would be dissipated 
during its passage to the thermometer well, 3. The pure 
liquid vapor generated in the percolator tube enters the 
boiling temperature well, 3j where a part of it is condensed 
along the thermowell and into the reservoir, 14* The remaining 
vapor goes to the condensation temperature well, 5, where 
again part of it is condensed and enters the reservoir through 
the condensate return tube, 15. The uncondensed vapors next 
pass to the condenser, 7^ and through the drop counter, 8, 
liquid finally returns to the reservoir through the return 
tubel2. This recycling proceeds until the equilibrium 
conditions - constancy of boiling and condensing temperatures- 
is attained. The differe.nce in the boiling and condensing 
temperatures is a very good test of the purity of sample under 
investigation. 

The thermometer wells 3 and 5 consist of 16 mm dia, 

13 cm, long and 13 mm dia, 14 cm. long, tubes respectively 
p 2 ?ovided with a glass spiral, 10, around it. The vapors ^ 

enter the top of the spiral, and travel to the bottom of the 
tube along the spiral thus providing intimate liquid-vapor 
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contact along the temperature well resulting in the measure- 
ment of "true” temperature. 

Temperature Measurement ; Temperatures are measured with a 
precision of 0,01 °C Toy use of copper-cons tantan thermocouples (3) 
and a Honeywell potentiometer. The thermocouples were cali- 
brated against a 25 ohm Leeds and Horthrup platinum resistance 
thermometer and K-3 potentiometer assembly. The calibration 
details are given in Appendix A, Each thermocouple is placed 
in the thermowell containing silicone fluid for better heat 
transmission. The four thermocouples for four temperature 
measurements, boiling and condensing temperatures of sample 
and reference substance, are connected to a key from which 
two main leads are connected to the potentiometer for the 
measurement of voltage difference between the hot and the cold 
junctions. There was some difficulty with this arrangement 
and hence later single thermocouple was used to measure each 
temperature ( as shovm in Eig,4). 

Pressure Measurement; The pressures are roughly measured 
in a long U- tube glass manometer, 10 mm i.d, and 15 mm o.d. 
containing triple distilled mercury. The manometer is approxi- 
mately 1500 mm tall which fixes the maximum operating pressure 
in the system. The correct operating pressure in the system 
was calculated at the boiling temperature of the reference 
substance (equilibrium temperature) from Reference (18), 

The error in pressure measurement is expected to be 0.1 mm or 
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less. ' ■ 

As o/iown in Fig.ij, each ehulliometer is connected 
to the manifold through a diffusion barrier-cold trap unit,, 
and also to tnree glass ballasts one large and two small,* 

The diffuoion barrier consists of a sintered disc. The large 
ballast consists of one 40 liter flask enclosed in a wooden 
ox to minimize the effect of outside temperature variation, 
connected to the system by a two-way glass stopcock, Si* 
Also, two small 4-liter flasks are similarly connected to 
the t\ro ebulliometers through two-way stopcocks, S2 and S3, as 
shown in hig.i, fijnction of the ballasts is to control 

and dampen the pressure fluctuations in the system. 

A Genco Hyvac vacuum pump is used for evacuating the 
system to 10 mm of Hg. Requisite pressure in the system 
IS maintained by bleeding nitrogen gas from a nitrogen cylinder 
through a calcium chloride trap. 

Each eballloEeter Is provided with sample entry tube 11, 
a stopcock, F;',wlth a ground glass joint, and a discharge 
tube with stopcock SIO.The sample could be charged Into the 

ebulllometor at S 11 without contaminating it with air and 
moisture. 


• « * • 
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CHATTER 3 
MATERIALS 

1 . Reference Substance Water ; 

Water was deionized, and distilled over potassium 
permang.anate and sodium hydroxide to give a conductivity 
less than 1 ,0 micro-mho , in the Metallurgical Engineering 
Department, IIT/Kanpur, The effect of any impurities on 
the measured boiling temperatures was within + 0.01 °0, 

Water was selected as the reference standard because it can 
be made in very pure form and accurate and well recognized 
vapor pressure data are available for water, 

2, Compounds Investigated ; 

n-Honanol, bromoform, 1-bromonaphthalene, and cyclo- 
heptane were AR grade samples. These were used as such without 
further purification. Benzyl benzoate, nonyl phenol, and 
propylene trimer were lab reagent grade samples and were 
purified by fractional distillation in a 20 plate sieve column 
of 2.54 cm, i.d., at a 5^1 reflux ratio. Each sample purity 
was determined by density and refractive index measurements. 

The densities were measured at 37.82Cin a pycnometer (22) and 
a Bausch and Bomb Abbe refractometer type 53-45-58 (Sr .Ho, CDS 64 6) 
was used to measure refractive index at 30°C, Table I gives 
the source, and n^® for each sample along with the 

literature values, 

A rough calculation of sample purity was done using 
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tlie follo^ving foraiula (ll) ¥/hich is based upon the elevation 
of boiling point in presence of an impurity. 








(4) 


where = mole fraction of the impurity 

= enthalpy of vaporization taken at bo, cal/mole 

2 

T-^ = boiling point of the impurity 
Tq = boiling point of pure compound 
R = Gas constant, 1.98717 cal. g. mole 

The value was taken as the difference between the 

boiling and condensing temperatures. To was calculated from 
the experimental vapor pressure data using the Antoine equation, 
value was calculated at T-j^ using the Antoine constants 
and the Ciapeyron equation assuming vapor phase ideality. In 
this manner the purity of each sample was determined and is 
reported in Table 1. 
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CHAPTEE 4 

EXPBBIMIAL PEOCEPUEE 

With reference to Pig. 1, the system was tested for 
leaks after proper greasing of the stopcocks, etc. Once 
this was done, the following experimental procedure was 
followed to obtain the vapor pres sure -boiling point data, 

1. The whole system was cleaned with acetone 
and distilled water and dried, 

2. S6, S9, S10, S11, S12 were closed and S1 , S2, S3, 
S4, S5 and S8 were opened. 

3. The system was evacuated for about 3-4 hours to 

mm Hg. (as tested by Mcleod gage). S8 was 
closed to again check leaks in the system, 

4. If no leak, S2, 33, S4 and S5 were closed, the 
sample and water were filled into the respective 
ebulliometers through S1 1 and S12 from the small 

rz 

saraple cylinders (50 cnr ) ’without exposing the 
samples to air and moisture. This was done in 
the following manner. The sample bottle stopcocks 
was closed and the bottle was attached to the 
ground-glass joint at S11, S11 was opened to 

evacuate the space between S1 1 and the sample 
bottle stopcock. The latter was slowly opened 
after starting the condenser water. The ebullio- 
meter was thus filled by the required amount of 
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the sample. This proced-ui?e was also followed 
for charging water to the ehulliometer , 

5. The lowest operating pressure was adjusted in 
the system hy slowly "bleeding dry N 2 gas 
through S7 8.nd S6. The lov/est operating pressure 
was fixed by the room temperature and the vapor 
pressure of water. S1 was kept open during an 
experiment for minimising pressure fluctuations. 

6. After the pressure adjustment was made slow 
heating was provided to the sample and water. 

This slow heating avoided any bumping. After 
some time, smooth boiling was observed. The 
boiling and condensing temperatures were noted 
every 10 minutes. It took on the average about 
5-4 hours to reach equilibrium i.e. the boiling 
md condensing temperatures remained constant. 

The eouilibriura pressure on the manometer, boiling 
and condensing temperatures of sample and water were 
noted . 

7. The pressure was increased to a suitable value 
approximately 50 mm higher each time and the 
whole procedure i.e. I'd. 5 and No, 5 above was 
repeated . 

8. In this fashion the entire pressure range was 
covered, the highest operating pressure being 
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about 1000 mm Hg, At about 1200 mm Hg. , stop- 
cocks started blowing out and hence higher 
pressures could not be studied, 

9. After completion of this work, sample and water were 
removed, the systesri cleaned and dried and the 
experiment was ready for another sample. 
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CHAPTER 5 
RESULTS 

The experiniental results are reported In Table II 
to VIII. Columns 2 and 4 are the average of the boiling 
and condensing v;ell millivolt readings. Columns 3 and 5 
are the temperatures obtained from the calibration results 
at each millivolt reading. Column 6 gives the system pressure 
calculated at v/ater boiling temperatures as given in Column 5. 
using the Antoine constants from Reference (18). The tempera- 
tures reported in column 3 are given to more significant figai’cs 
than the accuracy demands. The vapor pressure-boiling point 
data as given in columns 3 and 6 were fitted to the Antoine 
equation, Eqn.(l) and the Cox equation, Eqn.(2) by the least 
squares regression and optimization technique on IBM 7044 
computer. The flow chart and. the complete program is given 
in Appendix B. Columns 7,8,9 in Tables II to VIII give the 
deviations in temperature and pressure from the experimental 
data using the Antoine and the Cox equations. The Antoine 
and the Cox constants are also reported in the above tables 
along with the standard deviations of the fit. The values 

reported in the above tables at the normal boiling point 
were calculated using lqn.(3) as followss 

Antoine Equations ■ 
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EXPERIMMTAL VAPOR PBESSUEB-BOILIH& POUT DATA K)R 
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Oox Equation; 

= 2.302585 ET^, (6) 

2 

where log = a + bl, + cT, 

® cox b b 

Here is the normal boiling point in °K. 

‘I'he ibitoine and the Cox constants for all compounds are 
reported in Table IX. ' 
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GHAPTEB 6 

SISGUSSIOM ME COECLUSIOE 

The experimental method chosen in the present investi- 
gation has been demonstrated to give reliable vapor pressure 
measurements (12,25). Once the leahs were removed from the 
system, WG found it rather simple to take the data which are 
both precise and accurate, G are was taken where necessary to 
use high purity samples, an,d to introduce sample and water, 
into the system without exposure to air and moisture. One 
of the major problems faced in such studies is the procure- 
ment and availability of high purity samples. Our method of 
introducing the sample is not full-proof but .is found to 
be quite satisfactory. Overall, the temperatures are accurate 
to +0.01 degree but in some cases the accuracy is not as good, 
because of the decomposition of the sample as discussed later. 
Provided the reference sample is very pure, this method of 
pressure measurement is ver^'’ accurate. In the present investi- 
gation great care was taken to purify water and the errors in 
the temperature raea.surement due to impurities in water are 
expected to be less than 0.01°. This would give an error of 
less than 0,1 mm in pressure measurement, 

Por benzyl benzoate and nonyl phenol slight decomposition 
was observed as temperatures wore increased. Hence at higher 
temperatures, fresh samples were charged for each point. In 
this manner the decomposition was minimized. The final 
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estimatea of the a.ccuracy of the present data for each 
compound are given in Table X. 

Reliable data on these compounds are not available in 
the literature for comparison with our data. The comparison of 
older data with our values is shown in Table XI. In general 
the agreement around 760 mm Hg. pressure is satisfactory but 
the data deviate considerably at lower pressures. It is 
felt that the present data are closer to the "true" data 
than the older literature values. 
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TABLE X 


ESTIIIATED 

ACCURACIES OT 

' VA.P0R PRESSURE- 

-BOILING 


POINT 

DATA 


Compound 

A.ppro ximate 

Maximum Uncertainty 

Purity, io 

Temperature 

°C 

Pressure, 
mm Hs. 

Benzyl Benzoate 

99.55 

0.1 

0.6 

Bromofo rm 

99.96 

0.01 

0.1 

1-Bro mo naphthalene 99.88 

0.02 

0.1 

Cycloheptane 

99.97 

0.02 

0.4 

n-Bonanol . 

99.81 

0.01 

0.2 

Eonyl phenol 

99.83 

0.04 

0.4 

Propylene Trimer 

99.65 

0.07 

0.6, 


29 


TABLE XI 


COMPARISON 

OP PEESENT EXIERIMEI7TAL DATA 

WITH THE 


LITERATFEE 

VALLDSS 



Ben zyl 

Benzoate 


Pressure 
mm. Hg, 

Temperature, °C 

At, °C 

Present 

Ref. (15) 

59.4 

216.27 

220.0 

3.73. 

709.0 

319.46 

320.0 

0.54 

760.0 

323.25 

323.5 

0 . 25 , 

867.0 

330.65 

330.0 

-0.63 


1 -Bromonaphthalene 


* 


Ref, (15) 


138,0 

195.70 

210.01 

14.30 

182.0 

■ 206.74 

220.0 

13.26 

236.0 

217.91 

230.0 

12.09 

305..0 

229.79 

240.0 

10.21 

387.0 

241.66 

250-.0 


483.0 

253.53 

260.0 

6.4r 

601.0 

266.10 

270.0 

3.90 

760.0 

280.66 

281.1 

0.44 • 

911.0 

292.73 

290.0 

■ 
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TalDle Xi 

(continued) 




Bromoform 


PreSvSiire 
mm. He;. 

7’emperature. "=0 


Present 

Ref. (14) 

^t,°C 

100.0 

100.35" 

85.9 

-14.45 

200.0 

115.21 

106.1 

- 9.17, 

400.0 

00 

• ♦ 

127.9 

- 3,86 

760.0 

:j40.65 

150.5 

1 . 87 


Uonyl 

Alcohol 




Ref. (20) 


60.0 

128.8 

139.0 

10.13' 

100.0 

143.19 

151.3 

8.11 

200.0 

164.7-1 

170.5 

5.71 

400,0 

189.31 

192.1 

2.7'9 

760.0 

215.11 

213.5 

■ - 1 .61 
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A, B,C 
Sj B j c 

B. P., BP 

P 


E 

T 


t 

"■g 


A 


noimtclatite-f: 

Antoine Constants 

Cox constants 

Normal Boiling Point 

Heat of Vaporization,Kcal/g.niole 

"Vapor pressure, mm Hg,. 

Cas constant, R =1.93717 cal., K~^ g-mole“^ 
Pempcrature, °K 
Normal Boiling Point, 

Temperature, °C 

Molecular volume of gas phase ' 

Molecular voIueio of liquid phase 
Mole fraction of th© impurity 

Biffenence ( Exptl-Calc. ) 


• « 
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APMDIX.A 

T]iBMQ.GOlT?IjB QAIiIBRA'IIOH 

A Sargent Thermonttor Controlle'^ constant temperature 
■batli (+0*01°C) and Ethylene G-lycol (upto 130°C) and Maua Oil 
(upto 300°C) as hath fluids were used for the caiihration 
purposes* The circuit diagram is shown in Eigi34 The hot 
junctions of thermocouples TCI and ’TC2 wore dippev?. to a 
depth of 15 cm in the constant temperature hath fluid and 
tied close to the Platinum resistance thermometer. The cold 
junctions were placed in dewar fiasks containing a slurry of 
crushed ice and distilled water, Piie caution was taken to 
avoid the formation of air pockets in the slurry. The, 
calibration is given in Table HI, The measured data are 
correlated by means of the following equation* 

B(m7) = a + bt + ct^ t°C (7) 

The constant a, b,e v/ere calculated by P-osenbrock 
Search analysis on IM 7044 computer. The following ’’best" 
values wore obtained! 

b =3.7991025 x lOT^ 
c = 5.3526331 x 10“" 

The function orror of fitting = 1.98 x 10”'^ 

Eor actual temperature measurement in the system, Pig, 4, 
Eqn.(7) was used. 




Fi§. 3 - Thermocouple edUbration Circuit. 
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TABLE HI; CALIBMTIOI^ OE THERE/IO COUPLES 


Current; 2,99 mA 

D.C. ; 3 .Volts 

Reading 

Thermocouple 



Mo . 

Reading 

Current 



eini (mV) 

Radii 

Temp., °o 

1 

1.04775 

1.10282 

25.888 

2 

1.40645 

1.14157 

35.698 

3 

1.49155 

00 

— 

37.519 

4 

1 .63205 

1.16159 

40.7S1 

5 

1.86105 

1.18177 

45.908 

6 

2.04265 

1 .19746 

49.901 

7 

2.29455 

1.21979 

55.594 

8 

2.48055 

1.23720 

60.038 

9 

2.74175 

1.25883 

65.568 

10 

3.17955 

1.29759 

75.500 

11 

3,63245 

1.33609 

85.399 

12 

4.67935 

1.37168 

94.574 

13 


1.42092 

108.962 

14 

6.30655: 

1.53343 

138.726 

15 

7.32173 

1.60417 

157.583 

16 

9.11585 

1.72224 

189.304 

17 

10,1295 

1.78550 

206.433 

18 

12.1825 • 

1 .90708 

239.613 

19 

15,8535 , 

2.10639 

294.776 




APPEMBIX^B 


schematic DIAGrAm of the programme— ( IBm_7044, I it kAnpUp) 


********* ■ 

* start * 

^ * -Jf * -if- ^ ^ 

* ■' ■ 

"if 

##* ■if-if-if'if -Jf «Jf4f if-if-if *^"if ^-if ■Jf*******#^4f"if#*4f **"^-)f** 


* 

MAIN programme- 


* 

* 

input data. 


* 

*„ 

calculation oF 

T,p requires.. 

* 

* 

subroutines... 

TSjImPL AnO EmFTEM* 

* 


¥r 

' *■ 

.. ^ . . 

***** 

' •■ *** 

■ * ' . . . . . 

****************************************** 

* subroutine ANToIn * 

* calculation OF antoine constants # 

****************************************** 

... * .. 

* ■■ ■ 

* ^ ^ ■ 

.■ ****** ■ 

■ ' *** ■' 

. . * ' 

****************************************** 

* SUBROUTINE COX AND JANICE * 

* calculation of cox CONTaNTS ^ 

****************************************** 

■ ■ . ■ .* '.'. 

^ * ■ 

* 

***** 

*** 

*******************************'?^********** * ******** 

* subroutine rLO«^’<^WApT ) * *» * 

* , 'T;\ ■ v: ■ ."'•'..-A . -Jf-**#**** STOP 

* subroutine*, * ** '* 


>K 
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MAH PROGRAM 


MAII-i PRO (-'RAH' For calculcation of Vapor Pressiii’e-Seaperature . 
I'his program — Requires ouDroutinoB: iSAJIBL A EMRiEM . 


SUBROUi lRES 

AliTOIii Calculation of Antoine constemits. 

OOX — Calculation of error from given Antoine and Cox 
constants. 


EUITIOH 
JAl'IIG 
MAI I IT V 
MIRIMI 


VSCIOR 


Function program for Antoine and Cox equations. 

Calculation of error from given Cox constants. 

Matrix inversion by tlie Gauss-Jordan method . 

' Optimization method ( Rosenbroch Search ) 

for Calc. Antoine and Cox constants, 

— — Calculation of variations . 


AIIOIR AITI) GOZ Require Subroutine MAIIITV. 

I.!IiriT'IY -—Requires Subroutine VECIOR . 
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IISl OP NOIAITONS 


Symbols 

Descriptions 

Units or 
Equation 

ACS,BCS,CCS, 

C| vn 

ACW, BOW, CCW 

Constants of calibration 
of thermocouples for sample 
and water, respectively 

E=A+B1+C1^ 

AP 

Approximate pressure reading 
from manometer 

mm Hg. 

DTI, 1)12, DW 

Derivative of 1 w.r;t. E 
for samnle. 1 and 2, and 
water; Ireepeotively. 

1, °C 

BMPS, EMPW, 
ESMP1,2, and 

EW 

Data from potentiometer 
reading for sample, 1,2, 
and water 

mV 

N 

Number of sets for experi- 
mental data 


N COMP 

Number of compounds 


Mil, M12 

Letters in A-Pormat 


IIILE 

Name of compound 


1SAMP2 

lemperatures of samples 
calculated from equation; ■ 

E=ACS+BCSJ’ . 

CCS, 1^,1, 

1 

Condensing temperature 


2 

Boiling temperature 



« « • 










:»is 


mf: 




mm 






..'.a-;-' J • *.v 

•:"*' -M'' i 


VSisll’fv- 


READ ALPHA, BETA 
D^NUM OF VARl, 
»0 1*1 SL 


Cj) cj-i; 

\+ ^ id ; 

FIND X| SATISFYING 


NFL^+I 


l«eUE+l 


Cnfu.eq 
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■ * T ■ » , 

1,4,' ^ U.'** ^ ^ t ^ 

X I ui ^ V 

.C ' 

C HtL f""’ i\Co i\ M'‘‘<*. c T' v'k" 'CA'^C^^ a i I ', ur '^' A A w R ■ rA< E o o u r\ F “ i £ iM P E A 1 , 

C R E'L. A i i <« H IP , 

V' ■ i |■■'I^■^ J- i'\i P K lA' ■' :‘ti'i ^ I.M^. 'PR 0 '^R^i"!’v^^c>* • ■ ■ 

' 1 -v fA I'i i* u AA ! \i u.r i t" 1 u. r; ' 

■'C ', ,' 

G ' ' ■ ■ ■ ■ 

■ ■ i j'w] u f\j w X {j [\j ■ I X 1 L *“ f ^ Vv ( 3' i i'Mp 1 i -> w 5 G F 2 i -* A i 9 r ^ T- ^ / 

u i ,MU,MO X P,,s ‘-'1 I i I 9 ^ \ c ^ jO i u { ;;0 > 

u iM ,-'U, J f *1 / J 't ,/ 1 1 1'l / i f'^'i P 2 i 3 vJ ) f r"' X P I L i M } 

'v- Ol'U'i V,>U / J M l''iC ‘^Z (NJ ' . 

■ u<''!iP i)R M R. l^Ar / Lr'*iF'W ,( 3 ^ I'S'Em? o ^ 3 a ), ' , 

C. I i\ I u .J u i I w 1 r u RiM A 1 

, > M , 'J ( 'V ! 1 M ^ ^ i i i i. } ^ 

, 31 Pg;G-;m r i IX *131 1 11 in ) /IX* lii^- * 2X*2 0A6*7A5 IaI-k-/] X* i31 c ) ) 

32 i CrXlM ( 2-r n ^ 23 X 9 ro /\ p, p l A a *42 a » ’ 1 IR » a T u r roR RAFCQi- V 

IP :ifi«<G) :: *3X*lPJj/lXAl3l(lrP ] /2n a *; M'PlHb i^pF Mv 1;,' J^^L' 

4 A 1 ^ f — Error 1: ■■ Epi piV .£ ■' i '^aL., G, 2- - ■t-'^o.nROR 4^'* IX * Iip-h-. * a i-.i'^i.i' . r'- 

■'■3.,'TCAL C •-^T-Ei\R0R PCAL P,: 1! IGP ^P-llRz; n,lX*3r; - / 1 X'* 1 3 iX ) ) ^ 

5 3 OKM A 1 ( 4 * i 3 * I 9 • 3 * 1 1 2 • 5 * 1 2 • 3 * F 1 u • 6 > P 12 • 3 ^ F 1 1 * 3 * F 1 o • 6 * 1 X * ^2 * F 3 • A '* r 
110#3 *t 8*3 tr 12*3 ff 8 « 3 *A 2 } 

34 FORMATriX*131(lH^U/lX*13l(lHZ)/) 

'■ 33 P‘ ORMAI ( IX */ IX f x-FoR lrroP RcAOInG FrOiI POTEiRT i AL-MFtFR E(iAV3.*iiV 
A' '■ l'lX*-jf ' If't FHAp) 10 * ■>,*#£rR:GR '20 +-« 0 # 0005 */ 

2 1 X f 'X' . IF" E •GREAT ThaN EQUAL ■ 10 • •••ErP’OK I4, +« 0#0U5^^1. ■ 

36 Format (1X'*/71X*40TMnUArU DLvIATIoN of lRROP FOR**-5i*/lUX*^fTEMpERATUR 
■■ .* ■ 14 . For. 'SAMPLE : 1 *F12#3*^^ UEGrEc. C'^^/ 1 UX'*;*TLmpERATuRE For sample 2 

2 14*5*x JL.OP44 l*/iQX **^f'TLMpLRMl URc-' '..F q'F XAT^pCrEF) =**Fi2*5*'X 

,3 l^LuRuc V* X /// 3 X * Jriu AuQVi- T^'^t^uL. ivATA .TOOK, FR-OM FIOT EnJ 1 Al-^ME T 

, 4R[; ) 

G 

g .,' ,' input data' ' 

RLAD lifTCOMP 
1 F'ORMAIOIG! 

uO 11 ICu“i»NGUiMP 

RLAD 2 t NM > i T 1 1 L t U ) 5 I = 1 * NM ) 

4 F 4 '**' I'*! A T { I 0 f 1 2 A U / > 8 A 6 ) 

iRlmo j # ^ uLv * CC^ 

f\ L A D 3 > A C W * L> G w * C C W 

3 F ORMIAT { 3L15 ♦ o ) 

READ AfN 

4 FORMATS 13) 

READ 5 * < EW( I ) ^ I«l| H ) 

READ Df CEaMPK 1) f 1^1 fH) 
read 3*(ESMP2(3 )*I^ 1 *N) 

5 Forma t { sf i o # 6 ) ^ 1 

read ApU )'» 

6 F 0 R fM A T ' ( i 0 F 8 • 3 ) - . ^ i ^ 
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READ 7»PNT1»PNT2 
7 FORMAT (2A2) 

PRi NT OUTPUT 
print 50 

print 51. (TITlE( I) »I=I,NM) 
print 52 


►TArTInG calculate ESAmplE Tn TEmpErAtUrES 
SDTW=0. 

SDTSi=0= 

SDTS2 = Oss 
SDTWeO. 

DO 14 1=1 »N ■ 

ESM2=ESMP2( n 
EF = EWd) 

ESMl = ESMPl(n 

call TSAMPL(ACS,BCS,CCs,Ts 1 »DtS 1 ,E 5 Mi) 
call TSAmPL(ACS,BGS,CCs,TS 2 » 0 TS 2 »EsM 2 ) 
call EmETEmC ACW,BCW.CCw,EF,TW»P .DT w.Dp) 

EMFW( I )=EF 

EMFS( I )=ESM2 

PEXPTLCn=P 

TSAMP2(I)=TS2 

DT2(I)=DTS2 

DW( n=DTW 

SDTS1 = SDTS1+DT1( I ) 

SdTS|«SDTS2+DT2( I ) 

SDTW«5DTW+DW( I ) 

PRINT53.PNT1. I .AP{ I ) , EsMl , TSl ,D tS1 ,EsM2 , TS2 .DtS2 ,PNT2 .EF , tw ,Dtw,P » 


-Dp.PNTl 

continue 

PRINT 54 
s PRINT 55 

Istandard deviation 

AVG1=SDTS1/FL0AT (N) 
AVG2=SDTS2/FL0AT (N) 
AVGW=SDTW/FlOAT(n) 

$DMT1«0. 

SDMT2«0, 

SDMW=0. 

Do 16 I=1»N 

, SDmTi=SDmTi+(DTi ( I )-AVG l)**2 

I SDmT 2 »SDmT 2 + <DT 2 ( I )-AVG 2)»-!^2 

|6 SDmWsSDmW+.DWI I )-AVGW)«*2 
^ STDnTi=SQRT (SDmT1/FlOAT(n) ) 
STDnT2*SorT (SDmT2/F.L0AT<n1) 
STDnTW* sort (SDmW/FlOAT<N) ) 
print 56*StDNTl,STDNt^VSMHTW 
ANtOlN ' " ’ ' 

continue 

Wop 

END 



-o ‘ c, r ; O ■- w o. n n m 
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IBFlC TSAMF^.. 

calculation OF Tt-MPERA'fup 


?u. rRO?''' E;vih(f./v).,, P’^ TElT 1 0-ME TEp pEACt.mg 


bUupoUTiNt TMMPL(AC.bC.Ct.To,pTp,Ei„, 
UElIA-oorI ( tiC*-«.2-4*'«- [ AC-E3 m )-m-CC ) ' . 

TS=(-bC4-DELTA)/(2.^CC) 

DESM = CuEt:-03 " ^ ' ■ . ■ , ' 

II {EoM*OE*]_u* ) UL3M“U*aE i'p 

DTUsDESm/ (iiC+;>.')?CC-»-lS) 

return 

CMD 



K®V- 


i; 

v’’ 

■e 


^ ^ LmLILm Fop CAlCUlAT ^o^; pP VAnQp p 
WA f LR • . » . , , A*^loiNC EoUATION 

^UuROU 1 I Nu <,,Hl TEm( AC , uC , v-C jEF , TW >p ,UT W>Up ) 
DEL 1 A-iCOR T ( DC».#2-4. AC-FF ) 5J-rc) 
TW^M-BCt-DFLTA)/.(2.-!KC.J 


pESSUrE pELATlONSriJp 


IF^ C LI- .GL #iO* ) DESM = 0*5E-G2 

D T Vs' = U i-; f. : / ( B 'C,4' 2 • «-,C:C i'? T vV ) 


IV rr'.-;,(., i ..?8,3ao) cn in i 
AW -0*FlH42EAt h’ I 
U W ~ u • 1 '/ *.) I 'i u 4 
Cw-i^O « 2 "I A I UL U3 
GO To 9 

1 I F { r.w ,(jT * 886 i Go 1 0 2 
AWs=0»81:or-'.Ae9E 01 . 
BW-0*r/6/2A2E ,94. , ■ 

. CW=U. 236291: 0,3' ' : " ' ' 

■ , GO To 9 , 

2 .IF, F1v^*GT#4a.216) ,G0 To 3 
,,' 'AW=0:,,a0B86,T67.fTF0l "..'''..V ’ 

B:Wat)« 1 73939] E U 4 ' • F- 
-CWsO.p^Ki 0E,'03" ■; ' . 



46 




! W 


8 ‘ 

,9 


Cw=0.23410E O 3 
Go To 9 

IF f I'^»GT,*j9,<93g, j -qq . 

AW=0.80464202E qi , 

[iW=0.i, 715429E 64 
'-W = n,23214E 03 '' ■ ^ 

Go In q ' ' "" 

IF ^ TW*GT ,69,625 )' Go Tn 5 
AW = O,80n6295E 01 
UW=0,1695167E O 4 
CW»0,23O41E 03 
OC TO 9 ^ 

I F ( TW.GT ,70.920) Go Tn 6 
AW=0,79845588E 01 

oW*0,1678948E,-a4...., “'-.-vi-., 

CW«0#22897E 0$ 

GO To 9 . . . • a.v. ^ 

;I F { TW , GT * 89#7Q5 } ■ ‘45,' 
AW*O,79634280E : 

i^WsO, 16659246 

CWS0.22777E 03' /I ' , ,0;:, 

60 To 9 ^ ^ 

IF ( T W , G T * 99 ), |, 55 '....To^ 8 -'?l 4 S 

AW = 0. 79483960G-"'01' 

BW«0, 16563906 44 
CW»0,22686E''a3" ' ‘ii., 

GO To 9 ,■■■■' . ’ 1 'V ‘4|;- 

AW*0.79186968E'Oi . ..O*'"" I! , 

B W-0 . 1 636909€ ^ 

CW*0.22492G43 4' 
p.6Xp(2*SPl||f:*M" 
op»DTw#p#{i^r/(:aw,,*^ 

RETURN 

END . M 
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0),»ATCAl{3 

0* ♦A (10,10 

CoX.( 10). Aw 


'^■'"'’‘'Of'!/ JO 

CQ,'lM0ri/JA 
‘^CM'-'ON/AM 
nL'A'L L^fA: 








M - U 0 J « 1 iiM j 

jUHnuc ;'^ 

A ( 1 ♦,!< ) sA i K , I ) , : ., . ,: 

Con t I nue: a ;: 

;Do ' /:o K«;i tM C3MSi9S3$M: 

p<K. ' 

Uo J*i tipi ; 

C**0 {3*15 ' ' ’ . 


PKiNT iiO* S ; m- ;>,• / 
f'RiNT' 120 r~:i: 

r::S. 


CnCALdJ-^'Al'oij^ 


AH^iorte2*2a||M/i.r' 




w ■ » .... ^ 



wfMM 




^ ^ I J A V G ^ 1 4 ‘*‘- 2 

1 ( 1 ) -A VCT ) »» 2 

wTAmdard devIati 

( CMp/F(_0 j j j 

( .i»‘>T/P|^QAt ( Ml ) ) 


»STDfjT. 


)' 1 , ' ''■■ ■' 4 •„' ■, U‘Ci:;!v 

^ ANTOInC CfAM 5 TA'NTS 




en*vv 

‘Caul c-j'xc.i'lfok’^ 

5 L NO pa 2.. 10 A'#,! 

.i^tUNT 'iSftriNai 
'P|INT .147.':A^ 
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rc CuA 
%j u 

'■ 1 M 


r I NE cgX ( u ,sp »M J n 
M'-N^iON Y(30),F(10,3 u 


*3U) »A(lu,iO) »B( iO,i) 


C 


. ' ' ’ " ^ ^ « 1 1 w f ^ 0 \ i A r 

^uMMON/ JOHN/TE,vip (30 ) p^ )' 

'3fX,«L35= i^.EiS.S/STX.^b?- ! «»E15.0/57X, 

^*»Cl5.8/37X,-K.DiO= ''’‘"i5.8/57X,^fB8= #,E] 

2 2_ n nKif;f4p(|).}.273*l5 ' ■ . ■. , , ■. 

'’I I NO* OP CONOTanI'O 

i 1 * A -'i \ t i .... ■ ' . 


-M“ 3 »tMJ I 


W i t, ? i 


u 0 i 

4U 

j-im 


f''\ 


ih-J K.:i 


or ruNCTfON 


- - 1 ^ V « I 

‘ ‘ J ) «At.00iu { AloCiO ( P ( J ) / 7 ^o , >y a .-UA/'t (g J. j:,) 




( I-I) 


00 


it 0 


i 1-'-^ t J ( 1% # i J I’F I I I J 3 J J J 

M.tf t hit 


/ u 


i i t J } & i ^ 

J W I ^'4 ^ I (4 

H 1 *J3 -IN oTit-: 

j Q i* i 4 ‘J 1 I ?*'*! U |» 

Do 30 It, l,H 
Oq -jQ K~iti 
A ( is » { J a 0 , , , 

DO 0 0 J » 2 » N 

A 

Continue 

A i I » K ) «A { K * D ■' ■■• ', ■' ^ ■ . 

CCNTINUC ■; an , .■ 

D ‘0 7 0 2 »M ■■ •. A- 

U t j-0* 1 3 » j ■ ' ' .A ' ■ ., . ■„ 

OO 7 u .J*l *N 

P' i K » 1 5 * V ( j l »F. ( K* J ) -+0 { K f t I ' 

CONTINUE', ,, • ■' ' ■ '.Ai' 

*'• A L t. ?.', A I' I 'MV { A t, 0 f a; i h \ » ], 1 1 , j ' ', ; ''' 

7 I ,‘ i '1 , 1 2 0 » ( 0 S r », 1 ';! f ,I i-f 1 1 M J 1 8 P' : 


w.,a1.cD'L A I J ON 01 , ,,. 

1 1, U A,*} J 't C {, 0 » Op' ♦ <4;*' f }: 


10 


IND 


iNOE 
'fUPN ,.' 


AAyyWfeAAAlA 


iho~ 
>• S/ 



'54a; 
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k 


C ■ ' 

ov-'uKU'J jAHKL(U,Up#(i^,K} 

u i H'-N'-' i ON 1^*1 >»»■’£ 1U»3U) fDIFpOu) ■ / 

Common/ JO iiN/ I t;MP { 3u J , p { 30 j ■ ., ... V, . 

C0M:‘'A'iv/ JAhi£v)/N , '■ 

v.;;.,-';: ‘Ni/nLLi-N/T £ JO) ' ■ ’ , '"A '• 

I-' K I fi) ! 1 0 • . ■ 

i OKmm J { 1 X » 1 28 £ IH- ) /5x tKp ( ) MMiiwK-iipy ¥-Ti''mti\ rr r, - ' 

» JX»«0Ii P £ £-AL-00S) « .7X*4:-i OG/ A-'ny t «• «« * I r » 8X » ( GAL > MMHG 

U ^ ^ 

oo j«i . „ . G : ' 

wCAL“Ui» , ‘ ■... :. ' , - ‘ ■ . /': ' . , 

iJO A“i*M ■. ,■ ■ • G'"...,. 

wv.ALo-w'v.AL '*'U £ K # 1 ) jrF i K* jj . ■ ,, 

ACwA ‘i L Ai ' ( L • Jg S'OCAl-.) . ,■ 

Cb® £ 1 £ £ j } ) tACOX ■ '■ ■ '''V ' 

.PCAL» V6b«#bXP£2*i026#C0) - ■, ; 

0 1 1'-P < J 5 «PCAl-«P£ j ) ■■ ' ; 

•^1>P«;:>0P»-DI<-P£ J J - 

AlK'«l»/T£j3 - . ' ■ V, ■•■P’:. 


-PCAL'-PC J) 
■Olt'fUj} 

£ J 3 


80 


TUHTlJ}-2n0l5 ' 

f''(< i M r 1 A w # P < 4 ); * T- C 4 3 ; » pCAb;; p { j l § 

i '/■ F 1 8 » b 3 ' ' g' , ■ ;; ■ 

Cowr inul ■ ■ ; :VG:/r;:.' ./-p-G'/'-V, yv gp- 

,AVGP<*'i»Dp/FLOA'f ■' ' 


‘t'»A;CoX'»AT|i 


■OU (0 


-1*?^ 


7'& ■' bMP«b?4p-’i‘ £ w U'P £ 43 •*AV0p j''^fi'2’' ;f '.’ ,-',1 
r brop»04jKf (bHP/FL0AT|i^'')''|;: ■ ■■ ' 

PL .fFfN3‘ /l42*AVCjpi'bf&bL vp LP,''"'p ^p::p:..p vp'^'Pi "Pv' 

ORM*-* T' {./ )bK »'■*<' A V£R'A£a,«i, ,op -^v'l-F -p f'CAp.~o,^$‘ 3 

, ,iROR For RREObU.RL J ■ '■ '.V'.' 

;P'^'RU’URN ■ 
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I- i l,»f 1 w (•'lA i I V 


iu 

ij 

20 

JO 

40 
4 A 

00 
ua 
/ w 

{JO 

fj 0 

oU 

|U'.^ 

I, U 0 
ilo 

A ,Ao 
Ai+y 
i AO 
iA<tO 

1 /o 

(SOW 

2ofj 

210 

2^0 

2 .jo 

2AO 

2 f 0 

{J lo 

^20 
-;A .JO 

;;i *« m' 

.^'6^ 
■i' ? O' 

jfOw 

'420 

•♦40 

4&0 


i r 
OA 


u lil’n» ’PI''0T<20) .INDEx( 20,2) 

yc I thN = l r ‘•'0»>JR0K) . I leoLWj. jtQLUM I , ( AKAx ,T ,S'«AP I 

uu (io j*'l »i4 . .■■'■■, , ’ , ' 

i 'H v.j’M j j ko : ■ . . 

on A JO I --ifM ■ 
mKCm roli f’iVoT llLErtSfU' 

#^MV|aA«w#0 , , ■ 

y\j 4 05' i 

^ * io's »'60 ' 

00 ioO.. ■K*|,*'N ■ '■ ■ 'J ■ 

1 n +p i ooiA K,},-tA, :; J 

i i ( -Awii, { A'{ '■ 8 5 ♦loa'iioo'- ' ■ •' . ■ 

i AO»v» J , ■ ' .,. , ' . ■ ' •; ■' : , J" . ' \. A ■ 

iC 0 tO(M*iK , A ,. ' . •, A''A'A'Aa j' "’ A -'';. 

AJiAA-»A i J I A i A ' A' ■ A 

C<)Ml : ;A,' ' A'" , ■ ■' A ' ' ''■■ 

Coma iNUfc: " ■ , ■ ■ ^ .A' A‘:;, ' ■ ■ . ■, 

U>iV(VAt RA:)LOM;)'»iPlVoT(,ICO|.OM)+l '.'■•■'A' '' ' 

iwRfifCHA.No0,:«()40 ,Jo;.PU|. Pi V0|.. :E.t.'EMlN,ToA 0 NRIAGONAL' 

A .i4'O»260'»140 ' ' '' ' ' 

i ' A' a -A'. ^ -''A'A/^ 

WO 200 ' ;'R.*:;i':''i 4 :;‘A,,A' '-^A'A':'; A 

.W A i I ftOWtti'A - A. '.'../-s:; ■' '-A';' ■' A' ■ ' ' 

A 1 1 '■ "’a-'’’''; 

♦a iC0COM!RlA4i(j4p^.'AA4;-A/\;^ AA'^ ,A ".A; 

I i' ( j»i’l , 2'l»o>2&0''»£i4 '' aA'' '■'■■'A'' A^ ■■ 

00 (SSw 4«'l »4 a'^'A'A'A A'''''-:v,'4' :'A'4'"’' ''- A' ' '•';A''0A Va 

w iwAp «! w i' 1 ■«:#•£' );■■: "A ; "a ■■■ . / , 7 ■ : : : ' 

0 { i -04 *0 J »4,l i4o0OM4o£'4A a: A-C A 4'':'4 ''4 " 

»A‘ t RC^LP^^A»^■'<^^•o4A'pA ' 7 a'' A ' 7'7'.,.A':.A.'-,;' , .■■■; 

i NoiA j „;AA:j7A7'A^, ■'^'J'''' 7 a 4 a 7 :a 

1 ^4 O'E xvj i » 2 ■ )4i',co'£;uM '■ A ''' ■'7 ,:: '7 " ' :'7: a ; -'7 

i #Oii i|«A4iC0fcO'i^i:itd4£^ A-'^AAA;; :a7';a 

0 i, ’♦0 E ^ .p 1 .VO.T", 'fiiE'0:iNA‘ £Y 7;7p;l Vo-T; ■'££ Eptf wT,7' 'a' 7 

o i I £'PM » P'£7ii .pRM'tp;!''' V qI4 t V-^': 7,a' -a/'' .a-- 7:', ; '7 Av, ' 

4 l'Rj)LOM'plp0L4Wl4i'2:w77A'77''';7:7''';;' •,'A'.;'’''A:7 a7aaa,','''.a- 

M'Q '' A 'a' 77’' A-' '' •'■■•A" 

■ < f < ^ 2 A. 7 '^7a7a'';7A’4v'77A.,7aA'7?:A;v.:: 

'4'0 340 A£P|;7»Ma777'7c77a,v 77':a7a77A77.AAaA"7"'‘AaA''' 7 • ;!77a 

■ 4ii40ia44£7ii|i’i:l|0£yi:^p;M4W'Wt’';AA^ 

■: '■p44A‘£7£»>i7||7||A7':7'7'A7A77'7A''', A''7A;7AAA.AAA'''7;A-''’777':''r '7 A:77;A 

,7'f7'*4'l£lAii#7yf|7A|AAOAA:7':A7®^^^ 




"■'Siv'S' 

H , : 







''+6U uo 500 L«1,M ■' 

55A ^amJ 
■^^0 CONTINUE 

,,,,. cot 

^0 7lw 1,1,,^ ' 

03J ORO^^lNOEX-cinJ 

640 JCOlUm»IM01X{1 » 

J50‘00 705 

Iff i if 

Co,NTlNUr''-v-:‘,V" 

/io Continue 
75U Continue ‘ 

740 RETURN t''"" . 

: . End v”A,,,v,,lii|f:|ii 





THIS subroutine programme is For MINMI-ERROR function AnD 
parameters OE function model 

subroutine MINTmY{X»n) . „ 

OIMENSION A(10»10) »B(10»10) »P(10,10 »y(lO,lD) 

ItRdO) ,X(lQ).Y(lQ) 

C0MMON/JOHH/TEM( 3 O ) »pSV (30 ) 

: • COMMON/UAMES/NiC . . 

: DATA NS, MAxMIHrLlMlT, alpha, beta, E,FL/ 20d»a»3*0»0»5»(>» 


. NTRIALA-1 ; 

■■■■.■:N'EE.L.UR* 0 , 

'iTAAfT-),: ; „■ ■ . 

i:.''"'' 0-0, ; ; 

Ait>j:i*,p*o,;': 

I',,:;,, B( 

I;' "P'T IVJTA'QfP:;,' -'v'' 

h. ■ ■ .„''VT I , J'T'='P'<»'-. • ■ ' ' V’^ 

1 ?'^ IF rT 4 "Tsi;*v",: 


K.-l:' 

'■k 


:COfNTrNUE:'’;^:vdv V'"' 

■■F#FUnX'GnT G>N''T;" •' v; 
if(ntrial.lT,0) go To ,9 
IF(MAXMIN.E0.0J go To T5 
IFClIMIT.EQ.OJ go To 18 
IF(F.LE.Fm1 go to 13, 
IF(NSTA6E.lT.O) go To 5 

NFElUR»«NFELUR+1 , ■■ •, 

|FInGPXES.GT#OI Go 

E»»-^BETA*e 

NTRIAL*NTRIAL+1 ; • ' . 





SIMULATAnEOUS 
,D(10) , 0 ( 16 ) 

00005 ,0 */ • 









’y^ i . , ,,/* v,'v f^tTh 

*■> <- >' > 




PRlhit.' l'tP»o, 
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$IBFTC VECTOR 

subroutine vector( A,B,C,p,V,D,R,N) 

■ Dimension A(io*ioi»B(io, in) ♦C{io»id) 

do 1 J = 1*.N 
DO 1 1=1 »N 
DO 1 K=J»N 

1 MI »j)=A( I»j)+D(K)*Vn »k) 


Do 2 I=i*N 
B(I»1)=A<I=1) 

2 R(l)»R(li+BCl»l)«*2 
R { 1 ) =SorT {R f 1 ) ) 

DO 3 1 = 1 »N 

3 V(i,ij*Bt I=i)/R(i) 

• ’ DO 8 ' ' 

DO '4 

DO 4 T=1#N 

4 C{JilC)=etj»|(0+A( I* J)*V( I*K) 
DO 6 L=1»MJ 

. : DO 6 T=I»N 

P ( 1 1 J ) = p ( I ♦ u ) +c ( J » L ) # V { I ♦ L > 

6 B ( I »Jl=A:Ti»U)-P^I M ) 

DO, 7 

7 R tu) =R f U I »u ) «'*2 

..'''■ ,;.R:(0:l4$QRt''.tR,t Jll '■ 

: '.'DO 8' 

.,8 :v:<.i*-,jTs'S''t'i,*,uj./R(.ji ■;; ; 

" ■■RETUR.N 


p(10. 


>V(10,1Q) *0(10) ,R{ 10) 



Sl8FTC^;FUNX0W. 

F UNCT i on fun XON' ( 0 » N ) ' ^ 

' ■ ■, U1F1E^|$14N- .OFIv').'»X MU'I ■ ■ ;, 
Common/:46H:N/T{:30,) 

■ :'COM'MON/iA,M£S;/'NK. ■' 

IF tN,E0*'3J, GO To 20 
DO 1 Isl.N 
1 X ( n=o( I ) 

: ■ AX1 = XM )Oj0,*1 . ; , ' 

AX2*-X'{2T«1.E^03 
AX3»X t 3) M*Er-07 ,, 

AX4*X (4)-n.E 02 


SUMaO* 

'C'^:^:-"'DEFr:MAT4DN 'EUnCTIOM,,^; ■' 

./a;:'; I^UNK.. 

./:';/-''--.,r'^'-::'AuGP'te-ALOGlO (P { I ) /T 60 -. 








au ic roRQK’ . ^ , 

: Trii- t aNw PACCq,A’'':, CAt.CULATi; 

.. ■ iuuRo.M'fJric. Fa^KoR (x vm»f ■ 

'■; ' X U^as-PIFp ( 3D 1 »DiFJ 

■ CCHMC-*a JOHri/TX'sP.) sp ( 30 S 
.CCMMOM/JAi'^ED/MK ■ . , 

»\PRIKT output ' ■ ' ■ ■ . 

F0PPAT(2r20.3»2F2C*6 > 

WnTf fCRHATdXt^ PA^-^METiPO Dtl5..8,,* 
12 Format { iXf/ix,i 23 (MwiiO 
l»»&X»4pCAL MMHu ■' '■ 

'■ i'4 FoRMAThX»-:; BU.N'CTION ERROR =“»£l6 


pF ■-:?-»POR' oF r Af'iD T 


'0* («-s <■' 
w 

- 

5 

*5 nr 1 
- w 





D 

V 

'"I LC- 

ri~ ),'/ 

r-X 


-sp 

(MM 

iv% 

} 


0 

Ds ■" 


X? 


DI.F' 


mQ 

0 

- 

c 

-"\.L y ■” 

ON E 

FR 


n «. 

V 

f c 

1.6 

f 

8 

f 

3*' 

nry T 

?.T 

t 

n f'l 

It 

•')l' i 

2r 

1, 

,6 

#'8 


AVCpAGt D l,rp { o3S-^C>vL ) 
, , ; ': F;- ;; 

«•,.» # : ■ ,B f » F, 1' , 2 » * C: r? 


.e ‘'2i Format (iX>4 ANTn'l^iE. Cp^sOTAMTU *'• ,,=f , Fl’P r?^AtT',l; 

■ “• I5)v, '" . ' ■ pP A ' '■ '■ "■ ■■ , ' ^ A' A^;.> ' 

*25 ; Format c ix i k FOfiCT.iON. "eppor »ei6 *F a'-^-etd ■■o;£..v..a;'P":.p’; rFA;2;*:Fv«';'Mnw 
• uTu UEv'» T K'^;- »i- 12 • UE-G* ■■ ■, /,,,, A /‘p, V.- 


IF £ 

la. ED 


V 

GO"-T 

A'X-l- 

XOl) 


'0* 

I' - . . -■ 

A.X2fe 

-.XU’- 

i 

W't 
'' * 

*£-^03 

AX3 = 

XOD.)- 


»■ 

E'-:-07 , 

>< 

it 

X f 4 ) 


i-t 

£-■ 02" 

PR In 

T ■■ 13. 


'AX 



DO 1 . '■ . ■■ ■'''•■■■ ApvV/'G'r 

!'XtT( I ),a2T3:*:1£p,/.A^;P; p's/pppP'-;'' ■ 

wxAXl + AX24^^K^AAX2'#'f}<!•Xr«-2■^' 'Q ' 
A-1C',«-*B P;:,a'v 

C-;V“' (1 .~-A:XA/TkM' 

PCAL---760 , *EXP:{ 2 *,30:20:#C'T' '' /V/V,, vR:v:. '"' p/Ap 
U I FP ( IJ = p U ),'-P£Au’ '.'V, p:' ■ :F.; 

PR In'i' i3,»P i 1 5 ».T;i I ) »'-pCA,L;aDT;rp;T^T''T 

1 continufPi ' 'p^\^'"pp--''"Appv;pA;oi^pAA^^ 

SDrf'P~0*/ ; ■' p^p:;^Fp .,0 M.;p;:;p.Op';'Vp:M ^ 

.DO :3 I -1 i PFp ' ' ■ ■ -■ -.■■Fp.. 

3 ODIFppDPIFp-i^DlFPt.iy'' p 3 yAV ppAppR j 
'■ AypP--DPl FP/.FLO'AT.TNX 5 3:', 3 :F''F‘'':3’3';3'33p’ 

bMl''r:S«0* . A ■ . 33 ^^3.\p3V^3pp..'.■ ■'p>p3A' 

Do 6 . 1 « 1 » NK ' ' ' ^ ' ■' '3 ;’v A .:''^'p -'A'V; ^3^33 
0 ^mDpS-^OmDp.CjA {.', 0 i F p 0.1 T'O.AVbp.:J:'»#'2;33:,* 

0 TOV p = 0 Q R T '■ t OiMD.pS'AF'lq A.T ..OA lk.T: T'3''„:' pi-ipp;.: 

PR I MT . 1.4 » F »AV&P t OTD VM:Pp.:p-0.3,,v;|3V3Pp Pp^O 




,«w:i'3;f<' 1 « 'i(*'lsK®».M! t’ >F<lliS» 1 ) fi'«WC # , 


i's*- 




^ , V 





K , ./ • ■' 

■ Gn‘'flQ 26 ■'■-■., > 

B 'AlaXCD ' - ■ ' 

■ 02" ■■ ■■ 

r A'$*3^'{3 3«-l^E' -D,?. ■ . 

fe;. print 21»A1,A2*A1 

B: D,G- 22 1 = 1 >NK 

P . PCAL=EXP{2*3026*(A1»A2/(A?+T ( I ) 

p DlFpI i Tisp { D-^p.feAL 
I TCAl sA2/t AlMl;:CGlO(PU ) ) >-A3 
|22 DlFf{l)=TU)-T-AL 
I SDIFp^feO.- " ' , ■.■ 
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